Astronomy & Astrophysics manuscript no. atmospheres 'pub 


©ESO 2010 


June 25, 2010 





On the radiative equilibrium of irradiated planetary atmospheres 

Tristan Guillot 

Universite de Nice-Sophia Antipolis, Observatoire de la Cote d'Azur, CNRS UMR 6202, B.P. 4229, 06304 Nice Cedex 4, France 
e-mail: tristan.guillot@oca.eu 

Submitted to A&A, October 2, 2009; Accepted: June 4, 2010 

ABSTRACT 

Context. The evolution of stars and planets is mostly controlled by the properties of their atmosphere. This is particularly true in the 
case of exoplanets close to their stars, for which one has to account both for an (often intense) irradiation flux, and from an intrinsic 
flux responsible for the progressive loss of the inner planetary heat. 

Aims. The goals of the present work are to help understanding the coupling between radiative transfer and advection in exoplanetary 
atmospheres and to provide constraints on the temperatures of the deep atmospheres. This is crucial in assessing whether modifying 
assumed opacity sources and/or heat transport may explain the inflated sizes of a significant number of giant exoplanets found so far 
Methods. 1 use a simple analytical approach inspired by Eddington's approximation for stellar atmospheres to derive a relation 
between temperature and optical depth valid for plane-parallel static grey atmospheres which are both transporting an intrinsic heat 
flux and receiving an outer radiation flux. The model is parameterized as a function of mean visible and thermal opacities, respectively. 
Results. The model is shown to reproduce relatively well temperature profiles obtained from more sophisticated radiative transfer 
calculations of exoplanetary atmospheres. It naturally explains why a temperature inversion (stratosphere) appears when the opacity 
in the optical becomes significant compared to that in the infrared. I further show that the mean equivalent flux (proportional to T*) 
is conserved in the presence of horizontal advection on constant optical depth levels. This implies with these hypotheses that the 
deep atmospheric temperature used as outer boundary for the evolution models should be calculated from models pertaining to the 
entire planetary atmosphere, not from ones that are relevant to the day side or to the substellar point. In these conditions, present-day 
models yield deep temperatures that are ~ 1000 K too cold to explain the present size of planet HD 209458b. An tenfold increase in 
the infrared to visible opacity ratio would be required to slow the planetary cooling and contraction sufficiently to explain its size. 
However, the mean equivalent flux is not conserved anymore in the presence of opacity variations, or in the case of non-radiative 
vertical transport of energy: The presence of clouds on the night side or a downward transport of kinetic energy and its dissipation at 
deep levels would help making the deep atmosphere hotter and may explain the inflated sizes of giant exoplanets. 

Key words, extrasolar giant planets - planet formation 



1. Introduction 

Many decades before computers would bring the possibility to 
model in detail the evolution of stars, an analytical solution to 
the problem of radiative transfer applied to stellar atmospheres 
revolutionized the study of their structure and evolution. This so- 
lution due to the English physicist Sir Arthur Stanley Eddington 
states that in a static plane-parallel grey stellar atmosphere in lo- 
cal thermal equilibrium, there is a simple relation between the 
temperature of the atmosphe re and the optical depth of the radi- 
ation that it emits (Eddington 1916). 

Nowadays, detailed computer simulations are available both 
to find solutions for the complex problem of radiative transfer 
in stars and planets, and to study their evolutions. In particu- 
lar, since the discovery of exoplanets, models to predict or ac- 
count for their size have relied on more or less detailed radiative 
transfer models or tables (*Guillot et al.ll996HBodenheimer et al 



Fortnev et al 



20011: iBurrows et all [2003; Baraffe et al. 20^ 
2007 ", to cite just a few). Radiative transfer models in plan- 
etary atmospheres have highlighted the possibility of a bi- 
furcation of solutions, depending on the presence or absence 
of effic ient absorbers such as titanium oxyde and vanadium 
oxyde (iHubenv et aL _ 2003h. leading in some cases to tern - 
perature inversions (B urrows et alj|2007bl ; iFortnev et alj|2008l) . 
Observations of primary and secondary transits at different 



wavelengths have even brought the possibility to test these mod- 
els and led to the identification of key molecular species (e.g. 
Tinetti et al. 2007; Barman 2008; Madhusudhan & Seager 20oi 
iSwainetalJl2010n . 

However, the problem is complex and cannot be fully 
grasped by radiative transfer models that remain largely one- 
dimensional: this is because first the stellar irradiation field 
on the planet is intrinsically inhomogeneous and second at- 
mospheric dynamics plays a cru cial role in the glo b al en- 
ergy balance. In light of this, owman & GuillotI 120021) 
predicted that close-in giant planets should be character- 
ized by significant day-night photospheric temperature vari- 
ations, with the possilibity of an asymetry in the light 
curve due to heat transport by zonal winds. Both were ver- 
ified by observations in the infrared ([Harrington et al.l 120061 : 
iKnutson e t al. 2007a). Detailed calculations combining radia- 
tive transfer and atmospheric dynamics are now being performed 
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(e.g.'Cho et al.'2008ULangton & Laughlin'200 8l:[Showman et alJ 
2009; Menou & Rauscher 2009; Dobbs-Dixon et alj |2010t~to 
cite just a few), but are intrinsically limited by available com- 
puting power 

In parallel, it has been realized that a significant fraction of 
irradiated giant planets are too large compared to what stan- 
dard evolution models predict (e.g. Bodenheimer et al. 2001; 
Guillot & Showman 2002; Burrows et al. 2007a; Guillot 20ol 
Baraffe et al.u2008>;iMiller et alji2009) . The possibility to choose 
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tion of the stellar effective temperature T,, its radius and the 
star-planet distance D: 



Fig. 1. Planet receiving a flux oTj|^ from its parent star and emit- 
ting an intrinsic heat flux crT^^^. The labeled quantities corre- 
spond to radiative fluxes perpendicular to the atmospheric sur- 
face at the location considered. 0^ correspond to the angle be- 
tween the direction of incidence of the collimated irradiation flux 
and the local vertical, and fj.^ - cosfl,. Fluxes that are mostly 
characterized by visible wavelengths are drawn in blue. Fluxes 
in the infrared are drawn in red. 



as outer boundary conditions of evolution models between vari- 
ous atmospheric models, from those calculated with a maximal 
irradiation flux at the substellar point to those c alculated with 
an irradiation averaged ov er the entire planet (e.g. lBurrows et all 
500 3; Baraff^e et alj|2003h has left some confusion as to whether 
atmospheric properties may account for the discrepancy. 

These are strong motivations towards the derivation of a sim- 
ple analytical model to capture the impor tant physics of the prob- 
lem. This approach was akeady taken bv lHansenI (l2008l) who fo- 
cused his analysis on the observational consequences in terms of 
emission from the planet. In what follows, I will be mostly inter- 
ested in understanding how the atmospheric properties impact 
the thermal evolution of planets. I first derive the temperature- 
optical depth relation valid for an atmosphere which is heated 
both from below (intrinsic heat) and above (stellar irradiation) 
and compare it to some detailed atmospheric calculations from 
the literature. The study takes advantage of the fact that close- 
in gia nt planets sh ould generally possess an extended radiative 
zone jGuillot et al.iJ996.) . so that convection may be neglected, 
as a first step at least. I then study the consequences of the varia- 
tion of irradiation and advection on the mean temperature of the 
deep atmosphere. The resulting analytical temperature profile is 
used as boundary condition of evolution models, in an attempt 
to explain the inflated size of planet HD 209458b. Limitations to 
the models are detailed in section |5] 



2. An analytic radiative equilibrium model for 
irradiated atmospheres 

2.1. Setting 

The geometry of the problem is shown in fig.[r| The planet is re- 
ceiving an irradiation flux oTj^, cr being the Stefan-Boltzmann 
constant and Tii-r an effective temperature characterizing the ir- 
radiation intensity. In the cases to be considered, the incoming 
stellar irradiation can be considered as coming from a well- 
defined direction, with an angle 0» to the perpendicular of the 
atmosphere. (The angular diameter of the irradiating star as seen 
from the atmosphere is o' = arctan(2 * R»/D), being the stel- 
lar radius and D the star-planet distance. For an extremely close 
exoplanet D ~ SR, so that a ~ 22°, but in most cases of interest 
a < 10°.) To first order, the irradiation temperature is a func- 



(1) 



At the substellar point on the planet, the flux received by the at- 
mosphere is o-T*^_. On the day side of the planet, the flux perpen- 
dicular to the atmosphere is jJtCrT^^^, where s cos 6,. (Note 
that eq. ([T]i neglects the fact that the equator of the planet is 
slightly closer to the star than its poles, as this affects the irradi- 
ation temperature by a factor ~ Rp/2D, Rp being the planetary 
radius, i.e. generally by 1% or less.) It is useful to consider the 
irradiation flux avera ged over the entire pl anetary surface, that I 
note (tT^^, following ISaumon et al. l (fT996h . and define the equi- 
librium temperature as: 



- eq 
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(2) 



The maximal irradiation on the planet is crT*^ at the substellar 
point, the planet receives no flux on the night side, and the aver- 
age irradiation is crT^^, 1 /4th of the substellar point value. 

As indicated in fig. [T] at a given location in the atmosphere, 
a fraction AfitCrT*^ of the incoming flux is reflected, whereas 
(1 - A)fi,crT*^ is absorbed in the atmosphere and eventually 
reemitted back to space. A is the albedo, i.e. the fraction of 
the flux that is directly reflected. In general, this is a com- 
plex function of the properties of the atmosphere at the con- 
sidered location, wavelength and direction. While the incoming 
flux is relatively highly collimated, this is not the case of the 
reflected, reemitted and intrinsic radiation fields. Furthermore, 
both the intrinsic and reemitted fluxes are mostly characterized 
by long wavelengths (in the infrared), with characteristic equiva- 
lent blackbody temperatures that are generally less than 2500 K. 
On the other hand, the irradiation and reflected fluxes are mostly 
characterized by short (optical) wavelengths, with equivalent 
blackbody temperatures equal to the effective temperature of the 
parent star, ~ 5700 K for a solar-type star. As we seek an ap- 
proximate solution to the radiative transfer problem, this will be 
important because it shows that the thermal and visible radia- 
tions are mostly decoupled. 

I now consider a given location in the atmosphere and a 
ray of intensity /y^, v being its frequency and 6 its angle (or 
y u = co s 0) with r espect to the (local) vertical. Following lMihalasI 
(1978') (see also'Huben v et alJl2003"l) . I define three moments of 
the specific intensity as : 



(3) 



Jy is equivalent to the energy of the beam, AnHy is the radiation 
flux and 4n/cKy the radiation pressure. Given the parameters of 
the problem, 4-KHi„ - J 4nHydv - jJtCrT'^ on the day side (//, > 
0). 

The moments of the radiative transfer equation in a static, 
plane-parallel atmosphere in local thermodyna mic equilibrium 
and assuming isotrop ic scattering can be written (IChandrasekhaii 
ll960HMihalaslll978h : 



dHy 

— — = KyiJy - By) 

dm 
dKy 

— ^XvHy. 

dm 



(4) 
(5) 
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Furthermore assuming radiative equilibrium implies: 

y{Jy - By)dV = (6) 



/-•(XI 
I Ky 

Jo 



where m is the column mass (dm - pdz, p being the density and 
z the altitude), Ky is the true absorption coefficient, cry the scat- 
tering coefficient and;^^^ - Ky+a-y the total extinction coefficient. 
By is the Planck function. Using eq. (|4]i, the radiative equilibrium 
equation (|6]l can be rewritten as a conservation equation for the 
total flux, i.e.: 



H= I HydV^ 

Jo 



cte. 



(7) 



A closure relation is needed to solve the problem. lEddingtonl 
(119161) noticed that in the interior of stars, one could consider as 
a first approximation = A + Bji, implying J - A and K = A/3. 
Using KjJ - 1/3 in eq. (|5]l, and assuming a grey atmosphere 
Kv - Xv - Cte allowed solving the problem. 

Here I will follow the same approach but splitting the prob- 
lem into two parts: one for the incoming radiation, mostly in the 
visible, one for the outgoing radiation, mostly emitted in the in- 
frared. Strictly speaking, the approach is valid in the limit when 
the incoming radiation and the outgoing radiation are always 
well separated in their characteristic wavelengths. One may then 
solve separately the radiation field of both radiation sources, de- 
termine the source function, and hence the interior temperature 
profile. In practice, this is only partially true for heavily irradi- 
ated exoplanets: as the temperature at depth increases, the ther- 
mal emission is pushed towards the visible. In some case, the 
photosphere may be so warm that the contribution in the visible 
is not negligible. In other cases, the irradiation flux may be char- 
acterized by low eff'ective temperatures (e.g. if the parent star is 
an M-dwarf), so that the incoming and outgoing irradiation are 
not so diff'erent. However, in most cases we can consider that 
the region of the planetary atmosphere where the stellar irradia- 
tion is absorbed is characterized by temperatures that are signifi- 
cantly below the effective temperature of the star, so that the two 
radiation fields are mostly decoupled. This implies that charac- 
teristic mean opacities can be calculated for these two fields. 

2.2. The incoming (visible) radiation 

I first seek a solution of the radiative transfer problem for the 
incoming radiation in the visible. This is a long-standing prob- 
lem in planetary atmospheres, one for which the scattering of 
the incoming light by atmospheric particules is crucial in deter- 
mining the fr action of fl ux that is absorbed (e.g. Chandrasekhar 
I960- Meador& Weav er 1980; Toon et al. 1989; Goody & Yung 
19891) . For giant exoplanets orbiting close to solar-type stars (i.e. 
with orbital periods shorter than 10 days), the fraction of the irra- 
diation that is reflected back is generally very lo w, of order 20% 



or less, both from theoretical calculations (e.g. Sudars kv et al 



2003 



2008 



Hood et al..i2008.) and from observation s (e.g. Rowe et al 



It is interesting to note that since the visible radiation field is set 
by the stellar irradiation, can be calculated a priori. 

I now make an important simplification: I assume that the 
thermal emission from the atmosphere at visible wavelengths al- 
ways has a negligible contribution to the global energy budget 
and that one can hence assume that By ~ for v in the visible. 
The equations ^ and ^ can hence be simplified by integrating 
over visible wavelengths: 



dHy 
dm 

dKy 

dm 



(10) 

(11) 



This assumption is justified at low visible optical depth where 
clearly the incoming irradiation flux is intense and the contribu- 
tion from the atmosphere is comparatively small: Jy » By for 
V in the visible and t. <k 1 . It can be questioned deeper down 
in the atmosphere where most of the incoming flux has been ab- 
sorbed. This region is generally characterized by a thermal com- 
ponent of the radiation field that is much more intense than the 
visible part, so that for large enough t. values, both 7th s> J. 
and 7th s> J^^^-^^^^ Bydv. This justifies neglecting the visible part 
of the source function in this region as well. 

Following Eddington, I write fit = y/K./Jy. Note that this 
approach is valid in two extreme cases: for isotropic irradiation 
(e.g. Hu benv et al.112003 '). for which yu, = 1 / V3, or in the case 
of collimated visible irradiation (e.g. Meador & Weaver 198^, 
in which case //» - cosO,. Equations ( fTOt and ( fTTT i then write: 



dHjy,Hy) 

dm^ 



Xy 



Hy). 



(12) 



Because the incoming radiation is rapidly absorbed, it is fine to 
assume that both J. and H. vanish at great depths (m oo), 
therefore 



(Jy(m),Hy(m))^(Jy(0),Hyme 



(13) 



Again, it should be noted that 7v characterizes the intensity of 
the incoming radiation field, and that the thermal contribution 
from the atmosphere at visible wavelengths has been neglected. 
Furthermore, eq. ( flOb impUes that: 



Hy(0) = -//.7v(0). 



(14) 



In the case of an incoming radiation flux that is considered 
as fully isotropic, there is an inconsistency as we should have 
Hy{0) = -7v(0)/2. This inconsistency is at the heart of the 
Eddington approximation however It is due to the fact that the 
incoming irradiation flux cannot remain fully isotropic because 
of the larger absorption of grazing rays. Only an approximate so- 
lution can be found by neglecting the dependence on direction. 
In the collimated beam case, the solution is exact in the limit of 
no scattering. 



Snellen et al.ll2009l; lAlonso et a T] |2009l) . For simplicity, I 2.3. Tine outgoing (tliermal) radiation 



hence choose to neglect scattering: Xv ~ Kv 

I integrate the moments of the radiation field in the visible: 



iJy,Hy,Ky)= j (Jy, Hy, Ky)dV. 

«7 visible 

Similarly, I define a mean opacity: 



(8) 



Let us now consider the thermal part of the radiation field. As 
previously, we obtain average quantities by integration over ther- 
mal wavelengths: 



(B, 7th, ^th, ^th) 



Ky — 7y 



KyJydV. 



(9) ^th 



visible 



Jthe 



Jthe 



KyJydV. 



(By, Jy, Hy, Ky)dV. 



(15) 



(16) 
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As previously for the mean visible opacity, /fth is a function of 
temperature that can be calculated a priori for a known outgoing 
radiation field. 

The system of equations (|4]i to (|6l) is integrated over thermal 
wavelengths with the hypothesis that B ' thermal '^vBydv x Kx^: 



dm 
dK(h 



: A:th(^th - B) 

KthHth 



dm 

K±(J± - B) + Kyjy = 



(17) 

(18) 
(19) 



We can combine eqs. (fTST l. (fTTT i and ( fT9] l to find by integra- 
tion in m: 

//,h - //th(0) = -A<.7v(0) (l - e-^v".//..) _ (20) 

Separately, integrating eq. (fTSl l over column mass and using 
eq. (l20l l yields: 

^th - ^:th(0) = ^th [Hth(O) - fi,MO) (l - e-^^'"/"-)] dm. (21) 

Equations O and ([Till imply that //th(0) = // + /^.^(O), which 
allows integrating the above equation to find 

Kti, - K^(0) = HK^yn + MO)—^il (l - e-'^"''"-) . 



(22) 



Using the first Eddington coefficient for the thermal radiation 
field, fKth = KihjJih, we obtain the mean intensity: 



U = Ah(0) + H-f^ + m)-^ (1 - e-^^"^l^^) (23) 

jKth Ky jKth ^ 



2.4. Atmospheric temperature profile 

I now turn to the derivation of the atmospheric temperature pro- 
file. I first introduce a second Eddington coefficient to relate the 
thermal radiative flux to the mean thermal intensity at the outer 
boundary: ///th = ^^th(0)//th(0). I also define the optical depth 



T = Kt^m 

and, following iHanse nl (l2008h 



(24) 



(25) 



Using these relations, eqs. ( fT4l i. ( fT9l ) and ( l23l l. we can express 
the source function at each level: 



B ^ H 



1 ^ T 

fmh fxth 



fmh 7 fKth 



(- 








jfKthl 





(26) 



The fluxes H and //v(0) correspond to the imposed heat fluxes 
at the bottom and top of the atmosphere, respectively (note that 
//v(0) < 0, as it corresponds to an inward flux). 

In the case of the collimated incoming irradiation, H — 
crT^^JAn; Hy{Q) = -fi,o-T^^/47T. As previously for the visible 
flux, /jfth = 1/3 is valid both for a collimated and an isotropic 
radiation field. The second Eddington coefficient is generally 
chosen to be either fnth - 1/2 or fnth - 1/ V3- The former 
value is derived from the assumption of isotropy of the outgoing 
radiation field, whil e the latter can be shown to result fr om an 
isotropic scattering (IChandrasekhailll960l: lMihalasl[T978l) . The 




1000 1500 2000 2500 3000 
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Fig. 2. Temperature as a function of thermal optical depth ob- 
tained for different inclinations of the incident light, from (black 
lines, left to right) 0° (terminator) to 0, = 90° (subsolar point) 
[see text and eq. (|27]|1. Also shown are the isotropic average 
[eq. (|29l)1, 0° average [eq. ( l27T i1 and global average [eq. (|49] l1 (see 
text). The quantities used for the plot correspond approximately 
to the case of HD209458b, i.e. T^^ = 1469 K, = 100 K, 
A-th = lO^^cm^g"', ATv = 4 X lO^-'cm^g"'. The corresponding 
chord optical depth is indicated on the right axis [see eq. ( 1581 ) 
with H = 550km, R = 94370 km and 5=1]. 



1/2 value however yields a temperature that is closer to the ex- 
act solution at great depth. Using thus fKth = 1/3 and fnth = 1/2 
yields: 



3T\ 

int 



2 



3T^ 



2 ^i. 

- -H — + 

3 r 



(27) 



y__i^ 

7. 

The solution is equivalent to that obtained by iHansenl (l2008h . 
except for a different boundary closure relation, and the fact that 
eq. ( |27] ) accounts for atmospheric heating due to visible radiation 
(see Appendix). The t = limit for the temperature is 



T(T = 0) = <! ^-Tl 



, 1 T-4 



1 + 



7_ 



1/4 



(28) 



Thus, the temperature at the top of the atmosphere is lowest and 
equal to that of a non-irradiated atmosphere emitting a total flux 
cr(r^j + f-*^in:^ ^'^ly ^'^^^ when y — > 0. In all other cases, 
the absorption of part of the visible irradiation flux at high lev- 
els in the atmosphere pushes the temperature up there. In the 
limit of high values of y, the temperature-pressure gradient can 
even become negative, which is observed in models in the case 
of strong TiO/VO absorpti on (iHubenv et al.ll2003T:lFortnev et alj 
120081: iBurrows et al.l200 8^. (Note that however in this case, non- 
grey effects may dominate and should be considered). 

Wlien assuming isotropic irradiation, the incoming flux is 
written Hy(Q) = fa-T^^IAn, with / = 1 at the substellar point, 
/ = 1/2 for a day-side average and / — 1/4 for an aver - 
aging over the whole planetary surface (iBurrows et al. I I2OOI . 
Furthermore, we have seen in ^2.21 that in this case iz* = I/V3, 
therefore: 
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Fig. 3. Atmospheric temperature-pressure pro- 
file obtained for different irradiation angles 
(jj = cos 6,) and compared to calculations 
for HD209458b av ailable in the literature: 
lBarman~et iD ( l2005h (in green -who assume 
Tint = 500 K on the night side to Tj^t = 230 K 
at the substellar point), and IShowman et alj 
( l2008h (in red). The black lines are calculated 
as in fig.|2] except that Tin, = 500 K. The corre- 
sponding r = 2/3 optical depth corresponds to 
the P = 65 mbar pressure level. 



3r^ 



r 



V3 



y 

V3 



r 



V3j 



(29) 



Of course, the standard Eddington relation is recovered in the 
limit when Tirr <k Tint- For - /Cth, one gets that T'^ - 
3 r.'^^ /4(t-i- 2/3)-i-0.93 f ^im which is equivalent to the expression 
provided bv' Hubenv et alJ (l2003h in that case. 

Figure|2]provides a comparison of the temperature structures 
obtained for different values of the incident inclination 0, and in 
the isotropic approximation, assuming an incoming flux that is 
averaged over the entire planetary surface (/ - 1/4). The nu- 
merical values have been chosen as representative of planet HD 
209458b. Without advection and assuming a very small intrinsic 
heat flux, the temperature profile is found to vary dramatically 
between the substellar point (0, = 0°) and the night side of the 
planet. The temperature becomes mostly isothermal at levels for 
which stellar irradiation has been entirely absorbed and before 
the contribution of the intrinsic heat starts to become significant. 
At higher levels, horizontal temperature variations on the day 
side are smaller because grazing rays are absorbed efficiently. 
For grazing incidences (or equivalently near the terminator), 
eq. ( l27l i always predicts a temperature inversion. Alternatively, 
lowering the ratio /Cth/'<'v favors the formation of a temperature 
inversion even for vertical incidence. 

The comparison between the isotropic approximation and a 
0° approximation (i.e. assuming the atmosphere is at the substel- 
lar point but receives the irradiation flux that is the average for 
the entire planet) shows that the isotropic approximation yields 
deep temperatures that are smaller and temperatures at small op- 
tical depths that are larger. This is a direct consequence of the 
fact that grazing rays are absorbed at low optical depths. 

2.5. Comparison to models 

Dedicated temperature profiles at difi^erent incidences have been 
calculated for the atmosphere of HD209458b and are compared 
to the results of this work in fig. [3] The pressure was calcu- 
lated by assuming constant gravity g, yielding P - Tg/Kx^. 
Opacities in the visible and infrared were adjusted to obtain a 
goo d match to the more de tailed models of iFortnev et alj (l2008l) 
and lShowman et al at vertical incidence. As a result, this 



yielded Kth - lO^^gcm"^ and Ky - 4 x lO^-'gcm"^. The fig- 
ure shows that the match remains good for other incidences, and 
that differences between the analytic approximation and other 
models are of the same order as differences between the mod- 
els themselves. Note that the model by Barman et al. (2005) is 
a good match to the other solutions except at grazing incidences 
(fi - Q.l, 0.3) where the temperature is found to be much lower. 
One likely possibility is tha t the plane-pa rallel approximation 
used both in this work and bv lFortnev et al.l is overestimating the 
absorption at grazing incidences when compared t o the more re- 
alistic spherical-symmetry approximation used bv lBarman et al.l 

I now compare in fig. |4]the solutions obtained from eq. ( |27] ) 
in the cas e of vertical inciden ce (yU, = 1) to a detailed calcu- 
lation by iFortnev et alj (l2008l) for planets at semi-major axes 
between 0.025 and 0.055 AU from their star (assumed to be 
a solar twin). The planets have a 1 Mj mass and a 1.2 Rj ra- 
dius, a solar-composition atmosphere, and the irradiation flux is 
calculated as a mean on the day-side hemisphere only. Again, 
the values of the opacity coefficients were adjusted to obtain 
a fair match, which was obtained for k± - lO^^cm^g"' and 
Ky = 6 X 10"^ VTkrT^OOOKcm^ g"'. These coefficients are in- 
deed representative of values of the Planck or Rosseland means 
in these atmospheres, as is the weak temperature dependance on 
the opacity in the visible. They are compatible with the coef- 
ficients obtained specifically for the case of HD 209458b. The 
comparison also shows the limit of the model: at low pressures, 
the visible opacity is probably overestimated (except in the cases 
where TiO/VO are present), and at large pressures, the opacities 
are definitely underestimated. This is because collision-induced 
absorption and/or a rising electron abundance eventually pop in 
so that the mean opacities are to first order proportional to pres- 
sure (instead of being roughly independent of pressure). 

The relation with the same opacity coefficients can then be 
usefully compared to a solar-composition model for isolated 
planets/brown dwarfs bv lSaumon et al.l ([19961). In this case, the 
10 bar level is a useful value to tie the atmosphere (characterized 
by the part that is still at a relatively low optical depth) and the in- 
terior. This 10 bar level is also found to be within the isothermal 
layer for highly irradiated planets, so that it serves as a conve- 
nient outer boundary for the evolution models. Figure|5]provides 
the value of Tio as a function of Tint or Ti„, in different cases. 
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Fig. 4. Temperature profiles as a function of pressure obtained 
for different irradiation levels (corresponding to orbital distances 
between 0.025 and 0.055 AU) at vertical incidence (//, = 1), 
using eq. (IZTTi (plain, red lines) and comparison to the results 
obtained bv lFortnev et al.l (l2008l) (dashed lines). 
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Fig. 5. Relation between intrinsic or irradiation effective temper- 
ature and the temperature at the 10 bar level in planetary atmo- 
spheres. The line-styles indicate gravity: plain (g = lO-'cms"^) 
and dashed {g - lO'^cms"^). I assumed /cir = lO^^cm^g"' 
and Ky - 6 X lO^^cm^g"', except for the thin blue line in 
the Tint - 200 K case and variable 7^-, for which Ky - 6 x 



10"^(r/2000K)'''^cm^g-'. The results of the analytical model 
for Tin- - are also compared to those calculated for a solar- 
composition atmospheres of isolated planets (g = lO-'cms"^, 
diamonds) and brow n dwarfs {g = 10'*cms"^, triangles) by 
ISaumon et all (1 19961) . 



for the same values of /cir and a:v as previously. For the non- 
irra diated cases {Ti„ - ), one obtains a relatively fair match to 
the ISaumon et al ] (ll996h results. The model however separates 
from these numerical calculations both for temperature signif- 
icantly lower or significantly higher than the usual ~ lOOOK 
representative of giant planets at distances ~ 0. 1 AU and less to 
their star. This is due to changes in the mean opacities for these 
characteristic temperatures. 

Figure |5] also shows that in the case of a significant ir- 
radiation, the temperature is independent of gravity, and that 
7^10 Tiij. With eq. (IZTT i and fi, - I, this is easily explained 
by taking the limit Tint — > and t — » oo. In that case: 



' deep 



[4(7^ 3 



1/4 



(30) 



With the fiducial opacity coefficients for HD209458b y - 0.4, 
I find that the right hand side is equal to 1 .24, close to the em- 
pirical Tio/Ti„ = 1.25 used by Guillot (2008 i) on the b asis of 
models of irradiated planets calculated bv Jro et al.l (12005 ). This 
proportionality relation w as also shown to app ly to strongly ir- 
radiated atmospheres by iHubenv et al.l (l2003l) . who estimated 

rio/ri„. ^ 1.15. 



3. The non-uniform irradiation flux and the mean 
atmospheric temperature 

3. 1 . Consequences of advection 

I now turn to the problem of the non-uniform irradiation flux, 
and its consequence for the deep atmospheric tem perature that 
controls the global planetary evolution ( Guillot & ShowmanI 
l2002t lArras & Bildstenll2006t iFortnev et al.l l2007). As we have 
seen, the solution to the pure radiative transfer problem is a tem- 
perature field that is intrinsically inhomogeneous. For giant plan- 
ets in our Solar System, the combination of rapid rotation, long 



radiative timescales and of an intrinsic flux that is of the same or- 
der of magnitude as the incoming h eat flux yields a relatively ho- 
mogeneous temperature field (e.g. lIngersoU & Porcoll 19781) . On 
the contrary, close-in giant exoplanets should be locked in syn- 
chronous rotation (Guillot et al. 1996) and they are characterized 
by photospheric radiative timescales that are much shorter and 
intrinsic heat fluxes that a re up to 4 orders of rnagnit ude smaller 
than the irradiation fluxes jGuillot & Showmanl2002h . The ques- 
tion of the proper outer boundary condition to be used for evo- 
lution models is crucial. Planetary evolution models have been 
calculated either with boundary conditions inferred from calcu- 
lations applying to the whole atmosphere, to the day-side hemi- 
sphere or even to t he substellar point (see iBurrows et al]l2003l : 
iBa raffe et ai]|2003b . On the other hand, little attention has been 
paid to the consequences of the temperature inhomogeneities on 
the planetary cooling and of the validity of the different calcula- 
tions. 

In lGuillot & ShowmanI (l2002h . the problem of the planetary 
evolution with a non-uniform outer boundary condition had been 
approached by assuming that temperature differences may per- 
sist even deep down in the planet. In that case, an effective en- 
ergy transport from the day side to the night sid43 takes place 
simply to homogeneize the specific entropies at deep levels. As 
one would expect, the non-uniform outer boundary is found to 
lead to a faster cooling and contraction. 

However, th e time-dependent radiative transfer models by 
liro et a 0(12005') show that the radiative timescales increase very 
rapidly with depth, so that any remaining non-synchronous ro- 
tation or slow advection is susceptible to provide a very homo- 
geneous temperature structure at depth. In this work, I will as- 
sume that there is a level, deep enough in the atmosphere/interior 
at which the temperature is independent of latitude/longitude. 
(Note that this should be deeper, peharps considerably, than 



' Note that I use "day side" and "night side" for simplicity, but there 
is a strong variation in irradiation when moving from equator to pole 
that should be considered as well. 



6 



Giiillot: Irradiated planetary atmospheres 



the level at which the irradiation flux has been completely ab- 
sorbed). I hereafter turn to the derivation of the temperature pro- 
file in an atmosphere that advects heat horizontally. 

3.2. Radiative transfer solution witli advection 

I now consider that for each atmospheric location (6, (p) defined 
from the substellar point, mixing tables place by horizontal ad- 
vection and transports heat with a flux qWT. The radiative equi- 
librium equation becomes: 



K(Jy - By)dv = qVT 





or, 

/CthC-Zth -B) + KyJy 



qVT. 



(31) 



(32) 



The first moment of the radiative transfer equation (eq. |4]l 
becomes by integration 



Jo 



dHy 

-dv = qVT, 

dt7i 



and hence 
H(m) = H{oc) 



-f 



qVTdm. 



Note that since we envision that VT — > when m —. 
implies //(oo) = crT^^J(4n). 

Now, the equation for //th becomes: 

— = -Kyjy + qVl 

dm 

and by integration 

qVTdm'. 

i/th(0) = HiO) - Hy(Q) and therefore 

/-*oo 

i/th(0) = i/(oo) - I q'^Tdm' + yu.7v(0). 
Jo 

Inserting this relation into eq. (l36l l yields 
Htt = H(oo) + ^ij^(0)e-'^"'l''' - I qVTdm' 



(33) 

(34) 
oo, this 



(37) 



We now integrate the equation for the second moment of the 
radiation field: 



^th - ^th(O) = Hi<x,)Kim + MO)— 1^1 (l - e-^'"^"-) 

poo 

— \ Ki\i \ qVTdm'dm' 

Jo J;n' 

and by integrating by parts: 

^th = ^:th(0) + H{o.)K^m + 7v(0)— //J (l - e-'^"-/''-) 

Ky 

qVTdm' - k± I (m' — m)qVTdm' 

Q Jo 



(39) 



The relation for 7th can then be found simply from the first 
Eddington coefficient f^th - K±/J±. Then, using eq. (|32] | yields 



B = 7th(0)-H//(oo)-— H-y^(O)— — 1 - 

jK\h Ky jKth ^ 



—Kxhtn I qVTdm' - K^h I (m' - m)qVTdm' 
Jo Jo 



+Jy(0)—e-''^'"'^'' - qVT 

Kth 



(41) 



We use the relations fmh = //th(0)/Ah(0), HyiO) = -fi Jy{Q), 

/-*oo 

i/th(0) = //(oo) - Hy(0) - I q'^Tdm, (42) 
Jo 

and //(oo) = crT^^^, Hy{0) = jitCrT^^ to find an expression for the 
temperature profile at each location (t,/z, (p) in the atmosphere: 



34 



1 



3/Hth 3/A:th 
1 



3/Hth 3r/; 



7 

3yu, 



3r./ 



Hth 



-- |(— -H — ) r qVTdti 



n-- 

/ATth Jo \ «1 



l\qVTdm' - qVT 



(43) 



The relation is a complex one and its resolution goes beyond the 
scope of the present article. 



(35) 3.3. Mean atmosplieric temperature 

We are mostly interested in the deep atmospheric temperature. 
As discussed, in the presence of an efficient-enough advection 
process, the temperature at deep levels should become latitudi- 
/^gN nally and longitudinally homogeneous. I therefore average over 
latitudes and longitudes (defined from the substellar point) to 
obtain a global mean temperature that depends only on depth r: 



diud(p. 



(44) 



For a conservative advection scheme (in particular if q does not 
depend on yu, </> or T), j> q^TdcD - 0. (One could easily show that 
this would also be true of heat diffusion, as long as the heat flux 
(-33-) is conserved horizontally.) This leads to a great simplification of 
eq. ( l43T l which becomes after integration over all latitudes and 
longitudes (using yu, - fi): 



_ int 



1 _^ T 

/nth //fth 



1 

+ 

2 4 



Jo 



1 yU 


{-- 


' \ 




fnih yfKth 


I A' 


yfKth 1 





djs. 



(45) 



Note that we integrated the intrinsic flux over the entire planet, 
whereas the irradiation flux is of course integrated only over the 
dayside hemisphere. 

The integral term can be rewritten 



(40) Y = 



1 



1 



2/mh 3y/A: 



+ r 



(46) 
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or, in terms of exponential integrals E„(z) = t "e ~'dt, 

1 



1 1 

2///th 3r/i 



-E^ijT). 



(47) 



The E„ functions h ave 
jAbramowitz & Stegunlll964 : 



r//rth 

a recursive property 



1 



i«+l(z) = - [e 
n 



■zEniz)] 



which implies that with some algebra, Y can be written more 
explicitly: 



1 1 
+ 



2/Hth Sy/zfth 

+r(i- 



1 + 



EliyT) 



(48) 



With our choice of /j^th = 1/3 and fnib - 1/2, the equation 
for the mean temperature becomes 



— 3T. 



7"4 — 



2 

■ - + T" 

4 3 



3r; 



eq 



2 

3^ 



1 + 



2 

— h 

3 

+ y(l-y|£2(rr) 



(49) 



For yr » 1, the properties of the £2 function imply that 

2 



^^4 



37^4, 

int 



2 

- + T 

3 



-^-^ eq 



2 

3 3r 



This relation is very similar to what obtained in the isotropic 
irradiation case but with a (2/3y) instead of a ( 's/3)/y coefficient. 

This is the same equation as obtained for the day side when 
considering no advection, but replacing T^^ by T^^. There is 
hence a well defined mean temperature of the atmosphere at each 
level T that is independent of the advective process to transport 
heat from the day side to the night side, as long as this advec- 
tive process is conservative and takes place horizontally, over 
surfaces of constant optical depth t. 

Deep in the atmosphere (r » 1), the ra diative time sca le be- 
comes very long, almost proportional to (liro et al.l2005h . This 
means that any slow advection and/or any slightly asynchronous 
rotation would ensure an almost homogeneous mixing. Ideally, 
the location where this homogeneous mixing occurs should be 
used as an outer boundary condition for the interior and evolu- 
tion models. Under these assumptions, only solutions in which 
the deep atmospheric temperature has been estimated with an ir- 
radiation flux averaged over the entire planet are energetically 
consistent. 



3.4. Comparison to models 

Figure |6] compares temperature profiles obtained for 
HD 209458b and relevant to the planet as a whole by var- 
ious sources. First, as also shown in fig.|2] it can be noticed that 
the global average calculation described by eq. (|49T l yields very 
slightly higher temperatures at depth, but is otherwise extremely 
close to the isotropic approximation of eq. (l29b . The analytical 
solution (using the same values of the opacity coefficients as 
previously) is a good match to the more elaborate calculations 
by Fortney_et al. (2005) and Iro et al. (2005). These are however 
ab out 200 to 300 K coo ler than the temperature profile obtained 
bv lBarman et al ] d2005h . The figure also shows the envelope of 



solutions obtained in a dynamical circulation bv lShowman et al.) 
('2009^ which is helpful to show the range of variability of the 
temperature profile in one particular model including radiative 
transfer and dynamical circulation. 

The cross in fig.|6]corresponds to the atmospheric boundary 
condition used in calculations of the evolution of HD209458b by 
iGuillotI (|2008). These yield a fast contraction of the planet and 
a radius that is at leas t 10% smaller than the measurements (e.g. 
iKnutson et al.ll200"7al) . The dotted line indicates the temperature 
profile obtained for a 10 times smaller visible opacity: because 
of a more efficient penetration of the incoming stellar flux, the 
deep levels are hotter than in the standard case by ~ 1000 /T. 
This is approximately the amount that is needed to account for 
the observed radius without invoking extra heat sources (see next 
section). As can be seen from the comparison of published radia- 
tive transfer calculations, this is outside the presently measured 
range of temperatures that is obtained by detailed models. The 
possibility that the visible opacity may be lowered that much 
is unlikely because even in the absence of efficient absorbers, 
scattering will have a non-negligible contribution. An alterna- 
tive possibility is however to increase the infrared opacity by 
the addition of greenhouse gases at high altitudes while keeping 
the visible opacity to its nominal value, thus yielding a small y 
value. 

I now turn to calculations focused on characterizing the day- 
side hemisphere of HD209458b, motivated by secondary eclipse 
measurements. These calculations, in which the incoming stellar 
flux was averaged over the day-side hemisphere only are com- 
pared in fig. Uj Particu larly interesting are the calculations by 
iBurrows et air(l2007bl) and iFortnev et al.1 (5007) with TiO and 
VO clouds that have a temperature inversion and were found to 
be compatible with the colors measured with Spitzer during the 
secondary eclipse of the planet. On the contra ry, other temper- 
ature profiles with no inversion (e.g. the Fortn ev et alJ (l2007h 
model with no TiO and VO clouds) were found to be incompat- 
ible. When applied to a day-side average, the analytic approxi - 
mation is found to be a good match to the iFortnev et al.l (l2007h 
model with the same values of the opacities coefficients as pre- 
viously. The Burrows et al. (2007b) model is colder by at least 
300 K, and can be more or less approximated by the analytical 
model with an order of magnitude increase of the visible opacity 
-thus yielding a pronounced temperature inversion-. 

It should be noticed that the temperature inversion that ap- 
pears required by observations is yielding lower deep temper- 
atures. This is naturally explained by the fact that part of the 
incoming stellar flux is absorbed at greater levels and thus does 
not penetrate deep into the planet. One possibility to be investi- 
gated and that is out of the scope of the simple grey models pre- 
sented here would be for the presence of significant non-grey ab- 
sorbers: with opacities that are strongly wavelength-dependent, 
the energy in the center of absorption lines is absorbed high in 
the atmospheres, but the lower absorption in the wings allows for 
the possibility of a deeper penetration of energy at those wave- 
lengths. In any case, fig. Q highlights the fact that the present 
observational constraints relate to relatively high atmospheric 
levels, not to the deeper levels used to tie the atmosphere and 
the interior models. 



4. Atmospheric properties and the sizes of 
exoplanets 

I now reexamine the problem of the sizes of exoplanets in light 
of this atmospheric boundary condition with two parameters, Kth 
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Fig. 6. Globally averaged atmospheric 
temperature-pressure profiles for HD 209458b. 
The plain black line corresponds to the result of 
a global average over all values of yu, including 
the night side [eq. l|49H. The dashed line is 
a result of the isotropic approximation for 
an irradiation flux averaged over the entire 
planetary surface [eq. (I29b1. Th ese are com- 
pared to similar c a lculat ions by iFort nev et al. 
2005h. liro et al.l ilOOS). and Barm an et al. 



20051) . The shaded region indicates the en- 
velope of all possible t e mpera ture profiles 
found by IShowman etlTl ( l2009h . The cross 
corresp onds to the outer b oun dary condition 
used bv lGuillot et all ( [20M) and lGuillotI {200i) 
to calculate the evolution of the planet. The 
dotted line represents the globally averaged 
temperature-pressure profile obtained for a vis- 
ible opacity reduced by a factor 10 compared 
to our fiducial value, and is representative 
of the atmospheric conditions that would 
allow explaining the evolution of HD209458b 
with no extra source of heat (see text). The 
parameters used for the calculation are as in 
fig.|2l except I used T,„f = 300 K. 
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Fig. 7. Atmospheric temperature-pressure pro- 
files for HD 209458b averaged over the day- 
side hemisphere. The plain black line corre- 
sponds to a day- side average from eq. l |49l ), 
with fiducial values for the opacities (see fig. 1). 
Profiles resulting from visible opacities that 
are increased ten fold, decreased by 1/2 and 
1/4 are also indicated by dashed, dash-dotted 
and dotted lines, r espectively. Calculation s for 
the dayside from iBurrows et d] ( l2007bl) are 
shown in blue, with squares indicating the in- 
ferred photospheric depth of Spitzer secondary 
eclipse measurements in 3 IRAC bands be- 
twee n 4.5 and 8 m icrons. Similar calculations 
from lFortnev et al . (2008) are indicated in red: 
three calculations are shown, with two models 
of TiO clouds, and without TiO and VO ab- 
sorption, respectively. The envelope envelope 
of all possible tempe rature profiles found by 
IShowman et al.l ( l2009h is indicated as a shaded 
region. The parameters used for the calculation 
are as in fig.|6] 



and Ky. As before, HD 209458b is our proxy. First I rederive the 
difference in radius between the model radius, the photospheric 
radius and the transit radius . The problem has b e en discussed 
before (Hubbard et al."200r; 'B urrows et J] 120031: iGuillot et al.l 
[2006; Burrows et al. 2007a), but is calculated in the context of 
our analytical atmospheric temperature profile. I then compare 
the evolution of the transit radii for the different boundary con- 
ditions to the measured one. 



4. 1 . Vertical optical depth & photospheric radius 

Our approximation of a planar atmosphere is equivalent to as- 
suming that the pressure scale height in the atmosphere H = 
-dr/d In P is infinitely small compared to the planetary radius, 
i.e. H/r <Si 1. In the case of HD209458b, assuming a perfect 
gas, a mean molecular weight fi - 2.3 a mean temperature 



r = 1500K and gravity g = 980 cm s^^, H = 'RT/^g x 550 km, 
for a planetary radius R = 94370 km. Therefore H/r x 6x 10"^ 
which is very small compared to other sources of uncertainties 
and justifies the planar approximation. I will therefore consider 
that g is constant in the atmosphere. 

In what follows, I will use the following notation: X(r) will 
denote a quantity that is evaluated at level r but that is assumed 
constant in the atmosphere. The independent length variable in 
the atmosphere will be denoted z, z = corresponding to a radius 
level r as measured from the planet's center (see fig.[8]l. 

By definition of r 



dT - -pKthdz 



(50) 



In the context of a plane-parallel, hydrostatic atmosphere, this 
can be integrated to show that 



g(r) 



P, 



(51) 
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Fig. 8. Geometry of the problem for the calculation of the chord 
optical depth. 



or equivalently to relate optical depth and altitude; 

c/t dz 



(52) 



Note that // is a function of gif) (assumed constant in the at- 
mosphere), but also for T{z) which can vary significantly. The 
above equation can be integrated to yield 

z = -H(r)~e{r, z) ln[T(r + z)/T(r)] (53) 
where Oir, z) is a non-isothermal correction that is equal to 
1 



0(r,z) = 



\n[T(r + z)/T(r)] T( 



I rT(r+z) 
TT TiT')dT'. 
in Jrir) 



(54) 



In the limit of an isothermal atmosphere, 0(r, z) - 1 . In all cases, 
eq. (ISJb may be used to evaluate the height difference between a 
calculated level and e.g. the photospheric level t = 2/3. 

4.2. Chord optical depth & transit radius 

When measuring the size of an exoplanet from a primary transit, 
the level that is probed is higher than the photospheric level. It 
corresponds instead to the level at which optical rays that are 
grazing, a t the terminator, have an optical depth close to unity 
dHubbard et al.i 120011) . Neglecting refraction, we thus define a 
chord optical depth for this grazing incident radiation: 



Tch(v, r) = I pKyds 

As shown by fig.|8] (r + z)^ - + s^, hence 



Tch(v, r) = 2 I 
Jo 



Z + r 



(z2 + Irzyn 



dz. 



(55) 



(56) 



This equation may be simplified with our plane-parallel assump- 
tion (z/r <s: 1). We further use eqs. (|50] | and (l53T l to yield 



Tch(v, r) 



2r 



1/2 ^T(r) 



rx 



dr 



11/2 



(57) 



With a new change of variable Z - - ln[T(r + z)/T(r)], we get 

dZ. (58) 



Tch(v,r) = T(r) — 

Kth \Hir) 



(2.7rr \"^^^ r*' 
W)) Jo 



[ne(r,Z)z] 



1/2 



Using eq. ( |53]) . we now rewrite eq. ( l58b at level r + Az: 

xf 

-'o [7T0{r + Az,Z, 



1/2 



(59) 



The height difference between the photospheric level for which 
T(r) - 2/3 and the transit radius for which Tch(v, r + Az) - 2/3 
is: 



H{r)e(r,Az) Ini — 

Kth 



Inr T{r) 
H{r) T(r + Az) 



1/2 



f 

Jo 



r ~ ni/2 
[7T6(r + Az,Z)Z\ 



dZ] 



(60) 



In the limit of an isothermal atmosphere, d(r, Z) - I, the integral 
is equal to 1 (the erf function evaluated at infinity) and the ex- 
pression of the chord optica l depth and height diffe rence reduce 
to the relations proposed bv lBurrows et al. N2007ah . Thus in the 
isothermal case. 



Az = H(.)ln^(^ 



1/2 



In the more general case of a variable atmospheric temperature 
profile, eq. ( |60] | may be easily resolved by iterations. 



4.3. Thermal evolution and sizes of transiting exoplanets 

I now calculate how the transit radius of an exoplanet is af- 
fected by the outer boundary conditions, and specifically, us- 
ing eq. ( |49] |. by a choice of /Cth and k^. Figure |9] shows the 
result of the calculation applied to HD 209458b, assuming a 
low helium abundance Y = 0. 24, standard evolution models 
(lGuillot&More]|ll995l: iGuillotI [2008.) and with values of the 
opacities that vary so that y - /Cy/zcth ranges between 0.04 and 
0.4. The difference between the model radius (at 10 bars) and the 
transit radius in the visible is calculated using eqs. ( |53] ) and ( |60] |. 
With our fiducial values of the opacity coefficients, the mod- 
eled size falls short of the observed value by more than 10%, as 
obtained before jB odenheimer et al. 2001; Guillot & Showmai 
2002; Burrows et alj i2003c iBaraffe et all 120031; iGuilloU i20d" 
Miller et al. 2009). Accounting properly for the transit radius 
and not the photospheric radius has a relatively small effect in 
that case. Decreasing the value of y (either by increasing /fth or 
by decreasing Ky) does help in reducing the mismatch, but an or- 
der of magnitude increase is required in order to reproduce the 
observed value (even though this assumes a low helium abun- 
dance and no central core). 

Clearly, this order of magnitude change of the opacity ra- 
tio compared to the fiducial value is well outside the range of 
possibilities spanned by present-day models: As shown in fig.|6] 
the corresponding solution (dotted curve) is 700 to lOOOK hot- 
ter than elaborate radiative transfer models predict. This remains 
true even when considering only the day side average (fig.|7]). It 
therefore appears that alternative solutions involving additional 
sources of energy or non-radiative downward energy transport 
are required to explain the inflated sizes of exoplanets. 
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Fig. 9. Radius of HD 209458b as a function of age for different 
values of the thermal and visible opacities. The fiducial values 
are k*^ = lO^^cm^g"', - 6 x lO^-'cm-g"' (bottom curves). 
Alternative models are found by choosing /Cy = k^/S, k± - k*^/5, 
Ky - a:*/ 10, respectively. The dotted curves correspond to pho- 
tospheric radii (t - 2/3). The plain curves correspond to transit 
radii (Xchord = .2/3). T he measured age and radius of the planet 
dKnutson et al.ll2007bl) are indicated with their error bar. 



5. Limitations & complications 

5.1. General approximations 

It should be stressed that, as the Eddington relation, the rela- 
tions that were derived are only approximate because they as- 
sume a given dependence of the radiation field with direction (its 
quasi-isotropy). The discrepancy with the exact solutions s hould 
however be small (see lChandrasekhad 1 9601; iMihalasI 1 9781) com- 
pared to the other sources of uncertainties of the problem. 

More importantly, because a grey approximation is used, this 
means that we cannot study the consequences of wavelength- 
dependent absorption. As discussed previously, the presence of 
strong absorption lines will lead to an efficient absorption of the 
energy in the higher atmosphere while allowing a deeper pene- 
tration in the wings of these lines. This cannot be captured within 
the framework of the simple model presented here. 

Using a plane-parallel approximation implies that solutions 
at very low inclination angles are probably crude compared to 
the true spherical geometry. This generally tends to yield an 
overestimation of the magnitude of the temperature inversion for 
low inclinations. 

Finally, departures from local thermodynamic equilibrium 
may be significant at relatively high altitudes and are ignored in 
the present study. This should be relatively minor however when 
concerned with the temperature profile at relatively deep levels. 

5.2. Convection 

In this work, I have considered only purely radiative solutions, 
with however the possibility of horizontal advection (i.e. ad- 
vection along surfaces of constant pressure and therefore con- 
stant optical depth). The comparison with available models of 
HD 209458b has shown that it is indeed a good approxima- 
tion to present models of the atmosphere of this planet and in 
general this should remain true for all heavily irradiated plan- 
ets as th ese are bound to possess thick external radiative layers 
dGuillot e t al. 1996; Guillot & Showman 2002). Howeve r as no - 
ticed by iBurkert et al.. (,2005,) and .Dobbs-Dixon & Lid (l2008h . 



depending on the efficiency of horizontal advection, convection 
should be present at lower optical depths on the night side of 
these planets. In this case, the problem is modified: the pres- 
ence of a non-radiative vertical transport of heat invalidates the 
hypotheses made in § 13.31 This implies that at optical depths 
for which convection is present, the mean temperature profile 
should depart from that predicted by eq. (l49b . Because convec- 
tion would occur preferentially in regions with the lowest pho- 
tospheric temperatures, the effect would be a more efficient loss 
of the internal heat, or equivalently, a lower mean interior tem- 
perature. This would generally increase ffie radius problem. 

5.3. Advection & variable opacities 

Another complication is through the likely modifications of 
opacities and cloud co verage with tempera t ure an d irradiation 
level. As proposed by I Showman & GuiU"o3 (|2002|) . the signifi- 
cant temperature variations in the atmosphere coupled to the hor- 
izontal (and possibly vertical) transport should affect the chem- 
istry of the atmosphere. This is particularly true for condens- 
ing species, which could form clouds in colder region, settle 
to greater depffis and be present in lower-than-expected abun- 
dances on the day side. This in fact may explain the relative 
lack of Na observed in HD209458b (Showman & Guillot 200^ 
liro et alJl200? ). The consequence of a variable opacity field is 
ffiat, even if advection is purely horizontal (in terms of pres- 
sure), the averaging performed in § 13.31 becomes invalid, as 
in the case of convection. This is because we cannot consider 
that advection proceeds on constant r levels. Again, the effect 
has been investigated dBurkert et al]l2005l ; iDobbs-Dixon & Linl 
120081; iDobbs-Dixon et al.ll2010h . but in simulations in which the 
interior entropy was held constant, so that the consequence in 
terms of heat loss or internal temperatures has not been quan- 
tified. Qualitatively however, we can notice that opacities that 
increase with lower temperature/lower irradiation levels do tend 
to allow heat to penetrate efficiently into the planet near the sub- 
stellar point while suppressing its loss in low-irradiation regions. 
This would favor a slower cooling and would therefore tend 
to decrease the discrepancy between models and observations. 
However, we can see from fig. Qthat even if we would consider 
that one half of the surface of the atmosphere (the night side) is 
not participating in the cooling, we are still a factor ~ 4 short 
in terms of opacities to explain the observations. This problem 
should be investigated further but it is unlikely that it can explain 
the size discrepancy by itself. 

5.4. Non-conservative advection, kinetic energy transport 

Finally, it should be noted that advection is not necessarily 
conservative, that waves and shocks may occur, and that en- 
ergy may be transported as kinetic en ergy instead of heat. This 
is a complex problem (e.g. iGoodman 2009) and I only men- 
tion here that, only of order 1% of the incoming stellar en- 
ergy needs to be transfomed into kinetic energy and dissipated 
at deeper l evels to modify the evolution of close-in exoplanets 
dGuillot & Sho wman 2002; Guillot 2008). Again, this would po- 
tentially alter the modeled atmospheric temperatures. 



6. Conclusion 

An analytic semi-grey model to approximate the structure of 
plane-parallel irradiated planetary atmospheres was derived in 
ffie framework of the Eddington approximation [see eq. d27li1. 
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The model is parametrized by /Cth the mean opacity at thermal 
wavelengths, and /Cv the mean opacity at the wavelengths that 
characterize the incoming stellar irradiation. As in the usual grey 
approximation, these opacities are assumed constant, however 
the thermal and visible opacities may differ. The relation was 
shown to agree with more detailed calculations in the r ~ 1 re- 
gion, both as a function of the incidence angle and as a function 
of the mean irradiation level. The model qualitatively explains 
temperature inversions as resulting from a higher opacity in the 
optical than at thermal wavelength leading to a partial absorp- 
tion of the irradiation flux at high levels in the atmosphere. It ex- 
plains the proportionality relation between the deep atmospheric 
temperature (e.g. at a 10 bar pressure level) and the equilibrium 
temperature seen in detailed atmospheric calculations (small de- 
partures from this proportionality are due to variations of the 
infrared to visible opacities with temperature). 

The model was extended to include variable irradiation and 
a horizontal advection of heat. In the case of a purely horizon- 
tal (on constant optical depth surfaces) conservative advection, it 
was shown that the mean flux is conserved, so that a mean equa- 
tion for may be derived [eq. (|49]l1. Assuming that advection 
homogenizes deep levels because of the increase of the radiative 
cooling timescale, this relation should yield the proper bound- 
ary condition for internal structure and global evolution models. 
The temperature that is obtained is shown to be extremely close 
to the temperature obtained in a one-dimensional radiative trans- 
fer model assuming isotropy of the incoming irradiation and a 
mean flux o-T^^ that corresponds to an average over the entire 
planetary surface [eq. (|29]|1. 

A comparison of the results of the analytical model and 
of various available radiative transfer models for the transiting 
planet HD 209458b shows that the deep temperatures (at pres- 
sure below about 10 bars) that are obtained are generally about 
~ 1000 K too low to account for the observed size of the planet. 
Matching the observed and modeled radii requires a tenfold in- 
crease of the ratio of the infrared to the visible opacity in the 
atmosphere. This appears to be unlikely but the possibility mer- 
its to be investigated further given the ensemble of possibilities 
that remain in terms of atmospheric compositions and opacity 
sources. Alternatively, variations in the opacities (with higher 
thermal opacities in cold regions, possibly due to condensation) 
and kinetic energy transport are possible means to explain the 
size discrepancy by slowing the cooling of the planet. Progresses 
should be made by directly coupling radiative transfer calcula- 
tions to global circulation models. 



Burkert, A., Lin, D. N. C, Bodenheimer, P. H., Jones, C. A., & Yorke, H. W. 

2005,Ap/, 618, 512 
BuiTows, A., Budaj, J., & Hubeny, I. 2008, ApJ, 678, 1436 
Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007a, ApJ, 661, 502 
BuiTows, A., Hubeny, I., Budaj, J., Knutson, H. A., & Charbonneau, D. 2007b, 

ApJ, 668, L171 

Burrows, A., Sudarsky, D., & Hubbard, W. B. 2003, Ap/, 594, 545 

Chandrasekhai', S. 1960, Radiative transfer (New York: Dover) 

Cho, J., Menou, K., Hansen, B. M. S., & Seager, S. 2008, ApJ, 675, 817 

Dobbs-Dixon, I., Gumming, A., & Lin, D. N. C. 2010, ApJ, 710, 1395 

Dobbs-Dixon, I. & Lin, D. N. C. 2008, ApJ, 673, 513 

Eddington, A. S. 1916, MNRAS, 77, 16 

Fortney, J. J., Lodders, K., Marley, M. S., & Freedman, R. S. 2008, ApJ, 678, 
1419 

Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 

Fortney, J. J., Marley, M. S., Lodders, K., Saumon, D., & Freedman, R. 2005, 

ApJ, 627, L69 
Goodman, J. 2009, ApJ, 693, 1645 

Goody, R. M. Sl Yung, Y. L. 1989, Atmospheric radiation : theoretical basis 

(Oxford University Press, New York) 
Guillot, T. 2008, Physica Scripta Volume T, 130, 014023 

Guillot, T., Burrows, A., Hubbard, W. B., Lunine, J. I., & Saumon, D. 1996, ApJ, 
459, L35 

Guillot, T. & Morel, P 1995, A&AS, 109, 109 

Guillot, T., Santos, N. C., Pont, F, et al. 2006, Astron. irAstmphys, 453, L21 
Guillot, T. & Showman, A. P 2002, Astron. CrAstrophys, 385, 156 
Hansen, B. M. S. 2008, ApJS, 179, 484 

Harrington, J., Hansen, B. M., Luszcz, S. H., et al. 2006, Science, 314, 623 
Hood, B., Wood, K., Seager, S., & Collier Cameron, A. 2008, MNRAS, 389, 257 
Hubbard, W. B., Fortney, J. J., Lunine, J. I., et al. 2001, ApJ, 560, 413 
Hubeny, I., Bunws, A., & Sudai'sky, D. 2003, ApJ, 594, 1011 
Ingersoll, A. P & Porco, C. C. 1978, Icarus, 35, 27 
Iro, N., Bezard, B., & Guillot, T. 2005, Astron. & Astrophys, 436, 719 
Knutson, H. A., Charbonneau, D., Allen, L. E., et al. 2007a, Nature, 447, 183 
Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, T. M., & Gilliland, R. L. 

2007b, ApJ, 655, 564 
Langton, J. & LaughHn, G. 2008, ApJ, 674, 1 106 
Madhusudhan, N. & Seager, S. 2009, ApJ, 101, 24 

Meador, W. E. & Weaver, W. R. 1980, Journal of Atmospheric Sciences, 37, 630 
Menou, K. & Rauscher, E. 2009, ApJ, 700, 887 

Mihalas, D. 1978, Stellar atmospheres /2nd edition/ (San Francisco, 

W. H. Freeman and Co., 1978. 650 p.) 
Miller, N., Fortney, J. J., & Jackson, B. 2009, ApJ, 102, 1413 
Rowe, J. F, Matthews, J. M., Seager, S., et al. 2008, ApJ, 689, 1345 
Saumon, D., Hubbard, W. B., Burrows, A., et al. 1996, ApJ, 460, 993 
Showman, A. P, Cooper, C. S., Fortney, J. J., & Mailey, M. S. 2008, ApJ, 682, 

559 

Showman, A. P, Fortney, J. J., Lian, Y, et al. 2009, ApJ, 699, 564 
Showman, A. P & Guillot, T. 2002, Astron. & Astrophys, 385, 166 
Snellen, I. A. G., de Mooij, E. J. W., & Albrecht, S. 2009, Nature, 459, 543 
Sudarsky, D., Burrows, A., & Hubeny, I. 2003, ApJ, 588, 1121 
Swain, M. R., Deroo, P, Griffith, C. A., et al. 2010, Nature, 463, 637 
Tinetti, G., Vidal-Madjar, A., Liang, M., et al. 2007, Nature, 448, 169 
Toon, O. B., McKay, C. P, Ackerman, T. P, & Santhanam, K. 1989, /. Geophys. 
Re.v., 94, 16287 



Acknowledgements 

I thank J. Fortney and T. Barman for discussions on the charac- 
teristics of planetary atmospheres, the CNRS program Origine 
des Planetes et de la Vie and the Programme National de 
Planetologie for support. 

References 

Abramowitz, M. & Stegun, I. A. 1964, Handbook of Mathematical Functions 

with Formulas, Graphs, and Mathematical Tables (New York: Dover) 
Alonso, R., Guillot, T., Mazeh, T, et al. 2009, Astron. b- Astrophys, 501, L23 
ArTas, P & Bildsten, L. 2006, ApJ, 650, 394 

Baralfe, I., Chabrier, G., & Barman, T. 2008, Astron. Astrophys, 482, 315 
Baraflfe, I., Chabrier, G., Barman, T. S., Allard, R, & Hauschildt, R H. 2003, 

Astron. &- Astrophys, 402, 701 
Barman, T. S. 2008, ApJ, 676, L61 

Bai-man, T. S., Hauschildt, P H., & Allard, F 2005, ApJ, 632, 1 132 
Bodenheimer, P, Lin, D. N. C, & Mai'dhng, R. A. 200\, ApJ, 548, 466 



Appendix: An alternative derivation 

The temperature-optical depth relation described by eq. ( l27b was 
derived using two Eddington coefficients, fKth = ^th/-/th = 1/3 
and fnih = ^^th(0)//th(0) -1/2 equivalent to the assumption that 
the thermal flux remains isotropic, even at low optical depth. 
Another derivation is possible by imposing that the emergent 
flux should be equal to the sum of the intrinsic and irradiated 
fluxes (Mihalas 1978; Hansen 2008). In this case, one still as- 
sumes fKth = 1/3, but instead of using 7th(0) = //th(0)//Hth(0), 
we directly solve for yth(O) using eq. ( l23b and a relation obtained 
from integrating the equation of radiative transfer over thermal 
wavelengths: 

1 r°° 

/th(T = Q,n.,fi) = - B(t)e-'l^dt, (61) 
Jo 
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where /th(T = 0,fit,fi) is the intensity of thermal radiation emit- 
ted in direction jj from a point in the atmosphere which is irradi- 
ated by the star at an angle yu,. From eq. we obtain that 



f 

Jo 



[7th + yJv] dt. 



(62) = 



Using eq. (l23T l and integrating, one gets 
/th(0,//.,/i) = MO) + -^fi + 7v(0)-^ 

jKth jKthJ 
1 



+M0){7 



(63) 



We then impose that the flux emerging from the surface should 
be equal to the incoming flux, 4H + 4fi,Jy(Qy. 



f 

Jo 



(64) 



This allows expressing 7th(0) as a function of H and fitJy(O). 
Using then eqs. (|6]l and (l23T l. one gets after some calculations: 
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(65) 



This e quation is almost identical to the one derived by iHansenl 
However, a difference arises: the factor (3ii,/y - yln*) 
has replaced Hansen's {3fj.t/y). This is because the assumption 
of local thermal equilibrium implies that B = /th + yJv whereas 
iHansenl (|2008) assumes B - J±. Neglecting the yj^ term in 
the calculation of the source function implies that direct heat- 
ing from the irradiation flux is not considered: the atmosphere is 
heated only through the absorption of thermal radiation. In real- 
ity, the heating that is caused by the absorption of visible photons 
should be included, and it becomes a dominant source of heating 
when y/jd > V3. 

Equation ( |65] | is otherwise different in its form than eq. ( |27] ). 
but they have very similar properties. For example, at the sur- 
face, for an infinite penetration of the visible flux (y — > 0), then 
T'^{t = 0) — > (Tjl^j + jJtT^^)/2, ie the atmosphere still behaves as 
if it was transporting a flux (tT^^. = crT^ + crT* from below. 



As described in § 13.31 eq. 
the dependence on ju»: 



(l65l l may be averaged to remove 
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Solving this integral analytically requires more work than with 
the simpler temperature profile. The following relations are use- 
ful: 
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dx = -^ln(l +r)+ ^ln(l + - 
iy^^ 3 \ y 



3r \ 2y 




500 1000 1500 2000 2500 3000 
Temperature [K] 



Fig. 10. Temperature-optical depth profiles for different inclina- 
tions of the incident light. The figure is as fig. |2l but it com- 
pares the solutions obtained from eq. ( |27] ) (das hed black fines) , 
eq. ( l65T l (plain red lines), and the solutions of IHansenl (|2008|) 
(dotted blue lines). The thick curves in the middle correspond to 
global averages (see eqs. ( |49l ) and ( l67b respectively). 



The mean temperature profile is then shown to obey the follow- 
ing relation: 
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(67) 



Although the expression is more complex than eq. ( |49] ) which 
was derived using the second Eddington coefficient, the two are 
quantitatively very similar. 

A comparison of the various approaches is provided in 
fig. [To] Clearly, although the change in outer boundary condi- 
tion affects the form and complexity of the analytical solutions, 
the quantitative differences are extremely small, especially when 
compared to the large differences between published models for 
exoplanets (see § 13. 4b . Howeve r, there a re noticeable differences 
with the solutions provided bv lHanse B (I2OO8) for low ^l values. 
In particular, temperature inversions that should occur either for 
low yU or high y values are absent of Hansen's solutions, a direct 
consequence of neglecting the heating caused by absorption of 
visible radiation. 
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